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Abstract 

Planewave propagation in a simply moving, dielectric-magnetic medium that is 
isotropic in the co-moving reference frame, is classified into three different categories: 
positive-, negative-, and orthogonal-phase-velocity (PPV, NPV, and OPV). Calcula- 
tions from the perspective of an observer located in a non-co-moving reference frame 
show that, whether the nature of planewave propagation is PPV or NPV (or OPV in the 
case of nondissipative mediums) depends strongly upon the magnitude and direction of 
that observer's velocity relative to the medium. PPV propagation is characterized by 
a positive real wavenumber, NPV propagation by a negative real wavenumber. OPV 
propagation only occurs for nondissipative mediums, but weakly dissipative mediums 
can support nearly OPV propagation. 
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1 Introduction 

Although the complexity of the response of materials to electromagnetic radiation has long 
been recognized, the subdiscipline of complex-mediums electromagnetics (CME) came into 
prominence only during the 1990s — as Weiglhofer recently demonstrated in an extended 
review [1]. Besides nonlinearity, general complex mediums are not merely anisotropic, but 
are also magnetoelectric. Plasmas, ferrites, isotropic chiral materials, structurally chiral 
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materials, bianisotropic materials — all are excellent examples of linear complex mediums. 

Some of these materials occur in nature, others are artificially made. Add nonlinearity to 
the mix, and a bewilderingly huge palette of electromagnetic response properties emerges 
[2, 3]. 

While that stream of CME continues to flow unabated, a second stream of CME sprang 
in 2000. This new stream initially contained quite simple mediums: isotropic and dielectric- 
magnetic, but with the phase velocity of plane waves therein being in opposition to the 
time-averaged Poynting vector [4]. There has been some intermingling of the two streams 
[5, 6], but much of today's focus in the second stream is still on isotropic and dielectric- 
magnetic mediums. 

One current in that stream is the visuahzation of the negative-phase-velocity (NPV) 
scenario by inertial observers that are moving at a fixed velocity with respect to a certain 
medium (or vice versa!). This is a sensible issue, as human vehicles continue to penetrate the 
universe beyond our planetary atmosphere [7]. In a predecessor paper [8], we theoretically 
demonstrated that isotropic dielectric-magnetic mediums which do not support NPV propa- 
gation when viewed by a co-moving observer, can support NPV propagation when they are 
viewed in a reference frame which is translating at a fixed velocity of sufficiently high mag- 
nitude. Even more recently, we deduced the possiblity of orthogonal-phase-velocity (OPV) 
propagation — a phenomenon characterized by the orthogonality of the phase velocity and 
the time-averaged Poynting vector, when viewed by a non-co-moving inertial observer [9]. 

In this paper, we report a comprehensive study of NPV, OPV, and the conventional 
positive-phase-velocity (PPV) propagation of plane waves in a simply moving, homogeneous, 
dielectric-magnetic medium that is isotropic with respect to a co-moving observer. The plan 
of this paper is as follows: Section 2 is devoted to the Minkowski constitutive relations and 
planewave propagation. Section 3 contains numerical results for NPV, PPV, and OPV 
propagation in both dissipative and nondissipative mediums. Conclusions are presented in 
Section 4. 

A note on notation: Re {Q} and Im {Q} represent the real and imaginary parts, respec- 
tively, of a complex-valued Q. The complex conjugate is written as Q*. Vectors are identified 
by bold typeface and 3x3 dyadics are double underlined; v is a unit vector co-directional 
with v; the unit dyadic is /; and r denotes the spatial coordinate vector. The permittivity 
and permeability of free space (i.e., vacuum) are eo and //q, respectively; Cq = (eoyUj,)^^/^ is 
the speed of light in free space; cu is the angular frequency; and ko — cu/cq. An exp{—iu!t) 
time-dependence is implicit. 

2 Theory 

2.1 Minkowski constitutive relations 

Let us consider a homogeneous, isotropic, dielectric-magnetic medium, characterized by its 
relative permittivity and relative permeability iJ,r, in an inertial frame of reference E'. 
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That is, the frequency-domain constitutive relations 



D' = Coe^E' 



(1) 



describe the medium from the perspective of an observer located in E'. 

The inertial reference frame S' moves at constant velocity v = vv relative to another 
inertial reference frame S. From the perspective of an observer located in E, the medium is 
described by the Minkowski constitutive relations [10, 8] 



D = EoEra • E + 
mv X E 



mv X H 



B = 



where 



a — 



a 



Co 



q; / + (1 — a) vv , 



m 



1 — Er/IrP'^ ' 



1 — SrljLrP^ 

The transformation from (1) to (2) is achieved via 
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E' = (E • v) V + 

B' = (B . v) V + 

H' = (H-v)v + 

D' = (D • v) V + 
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[(I- vv) 'E + v X B] , 

. ^ vxE 

(|_vv) -B-— ^ 

[(|- vv) - H-v X D] , 

^ V X H 
(l-vv) -0 + —^ 



(2) 



(3) 
(4) 

(5) 



(6) 
(7) 
(8) 
(9) 



where P = v/cq. 

In the remainder of this paper, the propagation of plane waves in a medium described 
by the Minkowski constitutive relations (2) is investigated, with particular emphasis on the 
orientation of the phase velocity relative to the time-averaged Poynting vector. The reader is 
referred to standard works [10]-[12] for background details on planewave propagation within 
the context of special theory of relativity 
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2.2 Planewave propagation 

Plane waves in the chosen medium are described by the field phasors 

E = Eq exp {ik • r) 
H = Ho exp (ik • r) 



(10) 



Our attention is restricted to uniform plane waves, as given by the wavevector k = A; k where 
the unit vector k e M^. The wavenumber k is generally complex-valued, i.e., 

k = kR + iki , (11) 

where kji = Re {k} and kj = Im{A:}. 

To calculate the wavenumbers, solutions of the form (10) are sought to the frequency- 
domain Maxwell curl postulates combined with the Minkowski constitutive relations. Thus, 
we find [8, 10] 

k = ko - — 2" , (12) 

l-e(/?k.vj 

wherein 

e - '-f^, (13) 
A = l + (e,/i, - 1)5, (14) 

Note that the medium with constitutive relations (2) is unircfringent since the two wavenum- 
bers represented by (12) are not independent. The direction of planewave propagation de- 
termines the choice of sign for the square root term in (12). 

The amplitudes of the electromagnetic field phasors may expressed as [8] 

Eo = Ciei + C2e2, (16) 
Ho = 62 - a;eoer C2 ei , (17) 

where Ci and C2 are arbitrary constants. The orthogonal eigenvectors ei and 62 are given 

by 

ei = kxv, (18) 
62 = a X ei , (19) 

with 

^^^^ e(a,-k.v) ^ (20) 

c 

^0 



Hence, the time-averaged Poynting vector 



eip exp ( — 2A;/k • r 



( 




Pi + |C2|^cc;eoP2 



) 



(P) 



2 



(21) 



is delivered by (16)-(19), where 



Pi = Rc < — a 

[fir 

P2 = Re {e; a} 



} 



(22) 



The phase velocity 



— ° C- 



(23) 



is categorized in terms of its orientation relative to the time-averaged Poynting vector as 
follows: 

• positive phase velocity (PPV) is characterized by Vp • (P) > 0, 

• negative phase velocity (NPV) is characterized by Vp • (P) < 0, and 

• orthogonal phase velocity (OPV) is characterized by Vp • (P) = 0. 

As described elsewhere, the phenomenon of OPV propagation arises for nondissipative medi- 
ums when A < [9] . 

3 Numerical illustrations 

Let us now explore numerically the occurrence of PPV, NPV, and OPV propagation within 
the parameter space provided by {er,//r,v}. The velocity vector v is characterized by its 
relative magnitude (3 and its orientation angle 9 ~ cos'"-'^(k • v) relative to the direction of 
planewave propagation, in the remainder of this paper. 

3 . 1 Wavenumb er 

We begin by considering the wavcvector k. Notice from (12) that the wavenumber k is a 
function of the product e^/ir- In Figure 1, the real and imaginary parts of k are plotted 
for nondissipative mediums with e^/^r G (~10, 10), for (3 e {0.3,0.6,0.9} and the range of 
9 values given in Table 1. The imaginary part of k is null-valued across much — but not 
all — of the range €rHr > 0. 

Furthermore, at higher values of f5, kj ^ for larger positive values of e^/Xr- This reflects 
the fact that A can be either positive- or negative-valued for < e^/ir < 1, depending 
upon the magnitude and direction of v. Many schemes for materials satisfying the condition 
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< Eriir < 1 have been formulated [13, 14, 15], although we must note that dissipation then 
is hard to avoid. 

The real part of k in Figure 1 becomes unbounded for trf^r > 1, at points at which 
the denominator (1 — cos^ 6) in (12) vanishes. In view of our interest in the transition 
between PPV and NPV propagation, it is significant that changes sign as e^/^r increases 
from -10 to +10 only for 90° < ^ < 180°. 

For dissipative mediums, G C and fij. G C; hence, e,, = + ie^ and fi^ = + i/i^, 
where ef'^ G M and ix^'^ G M. The wavenumber is plotted function of G (-10, 10) 
in Figure 2 — for G {0.1|e^|, 2|e^|}, /^^ = =tl + O.li, (5 = 0.9 and the range of 9 values 
given in Table 1. Generally, kj is nonzero for all values of e^. As varies, strong resonances 
are observed in Figure 2 for both ku and kj. Notice in Figure 2 that the real part of the 
wavenumber changes sign as increases from —10 to +10 only for 90° < 6 < 180°. This 
behaviour is similar to that observed for the nondissipative scenario illustrated in Figure 1. 

3.2 Phase velocity 

Let us now turn to nature of the phase velocity, as delineated by Vp • (P). In Figure 3 the 
distribution of PPV, NPV, and OPV propagation is mapped out across the [36 plane for 
nondissipative mediums with e^/ir- G {1/3,2/3, 10/3}. The OPV regime arises for er/J^r < 1 
and for large values of (3, spanning a range of values of 9 which are symmetrically distributed 
about 9 — 90°. The NPV regime is restricted to ^ > 90° for both e^^^ < 1 and e^/^r > 1- We 
note that at e^yU^ = 1, the regions of NPV and OPV both vanish, in consonance with the 
Lorentz invariance of vacuum [11]. Planewavc propagation for the nondissipative scenario 
corresponding to e^/^r < (not shown in Figure 1) is OPV in nature for all values of /3 G (0, 1), 
since A < 0. 

For dissipative mediums, the distribution of the NPV and PPV regimes in the /39 plane 
is displayed in Figure 4, for ef^ G {±1/3, ±2/3, ±10/3}, G {0.1|e,^|, 2|e,^| }, and = 
±1 + O.li. When > and fi^ > 0, propagation is predominantly PPV in nature, with 
the NPV regime limited to 6* > 90° for large values of (3. Similarly, for < and fi^ < 0, 
propagation is predominantly of the NPV type, with PPV propagation arising only for 
9 < 90° at large values of /3. If and //^ have different signs, and the imaginary parts of 
Sr and Hr are relatively small, then the distribution of NPV and PPV propagation is evenly 
divided in the f39 plane along the line 9 = 90°. This even split between the NPV and PPV 
regimes breaks down if the imaginary parts of and/or /i^ become sufficiently large. 

We note, in particular, from Figure 4 that OPV propagation does not arise for dissipative 
mediums. Furthermore, there are no regions in the (39 plane for which the sign of Vp • (P) 
is indeterminate; i.e., both Vp • pi and Vp • p2 always have the same sign. 

In order to explore the orientation of the vectors pi and p2 relative to Vp, let us introduce 
the angles 



In Figure 5, contour plots of the angles 01 and 02 are shown in the P9 plane for = ±2/3, 





(24) 
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el e {0.1|ef I, 2|ef I}, and //^ = ±1 + O.li. We see that (pi < 90° and 02 < 90° for all values 
of j3 and 9. Thus, the projection of the time-averaged Poynting vector onto k is always 
positive. Therefore, NPV propagation arises when becomes negative-valued. 

The vectors pi and p2 are close to being perpendicular to k for large values of (3 when 9 
takes values close to 90°, in Figure 5. This is most noticeable when the imaginary parts of 
Er and are relatively small. We infer that nearly OPV propagation is achieved for weakly 
dissipative mediums for large values of f3 when 9 is close to 90°. It is clear from Figure 5 
that the orientations of pi and p2 with respect to k are generally very similar. In particular, 
for weakly dissipative mediums, the orientations of pi and p2 relative to k are almost the 
same across the entire P9 plane. 

The occurrences of NPV and PPV propagation are further explored in Figure 6, wherein 
the respective regimes arc delineated as functions G (—10, 10) and G (0, 10) for 9 G 
{30°, 150°}, 13 G {0.1,0.3,0.5}, and /x^ = ±1 + O.li. When 9 > 90° and (3 is small, the 
predominant type of propagation is PPV for values of much less than zero but NPV for 
values of much greater than zero. As (3 increases, the NPV regime expands, to the extent 
that when (3 — 0.5 most of the e^e^ plane supports NPV propagation. The reverse situation 
is observed when 9 < 90°: planewave propagation becomes increasingly PPV in nature as (3 
increases. 

The relative orientation angles (f)e of p£, {£ = 1,2), corresponding to the NPV and PPV 
distributions presented in Figure 6, are provided in Figure 7. Both (pi and 02 are acute 
across the entire e^e^ plane. Thus, our observation from Figure 5 that the projection of the 
time-averaged Poynting vector onto k is always positive, is further confirmed in Figure 7. 
Unlike the situation in Figure 5, the two vectors pi and p2 clearly have quite different 
orientations when viewed as functions of and e^.. Since f3 < 0.5 and 9 is not close to 90°, 
the propagation represented in Figure 7 does not approximate to OPV in nature. 

4 Conclusions 

From our numerical studies, for a simply moving, dielectric-magnetic medium that is isotropic 
in the co-moving reference frame, we conclude the following. 

(i) Whether the nature of planewave propagation is PPV or NPV (or OPV in the case 
of nondissipative mediums) depends strongly upon the inertial frame of reference that 
the observer is situated in. 

(ii) The projection of the time-averaged Poynting vector onto the direction of planewave 
propagation, as given by (P) • k, is non-negative. PPV propagation is characterized 
by ku > 0; NPV propagation arises when ku becomes negative-valued. 

(iii) OPV propagation only occurs for nondissipative mediums; for such mediums it occurs 
when the imaginary part of the wavenumber is nonzero. Weakly dissipative mediums 
support nearly OPV propagation for large values of /3, when 9 is close to 90°. 
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(iv) The orientations of the component vectors pi and p2 of the time-averaged Poynting 
vector are almost the same when viewed as functions of f3 and 6, but the orientation of 
Pi relative to P2 is sensitively dependent upon the constitutive parameters and Hr- 
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Table 1. Key for the values of 9 used in Figures 1 and 2. 
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Figure 1: The real (left) and imaginary (right) parts of the wavenumber k for nondissipative 
mediums, normalized with respect to /cq, plotted against e^/^r foi' G {0.3,0.6,0.9}. The 
values of d are identified in Table 1. 
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Figure 2: The real (left) and imaginary (right) parts of the wavenumbers k for dissipative 
mediums, normalized with respect to ko, plotted against for e {0.1|e^|, 2|e^|}, = 
±1 + O.li, and P — 0.9. The values of 9 are identified in Table 1. 
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Figure 3: The distribution of positive phase velocity (PPV), negative phase velocity (NPV), 
and orthogonal phase velocity (OPV), in relation to /3 G (0,1) and 6 G (0°, 180°) when 
erfir = 1/3, 2/3, and 10/3. The labels 'NPV and 'PPV shown hold for > and /x^ > 0, 
but must be reversed for < and /x^ < 0. 
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Figure 4: The distribution of positive phase velocity (PPV) and negative phase velocity 
(NPV) in relation to /? G (0,1) and 6 G (0°, 180°), for G {±1/3, ±2/3, ±10/3}, el G 
{lO.lef |,2|e^|}, and //^ = ±1 ± O.li 
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Figure 4: Continued. 
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Figure 6: The distribution of positive phase velocity (PPV) and negative phase velocity 
(NPV) in relation to G (-10, 10) and ef. G (0, 10), for 6 e {30°, 150°}, (3 e {0.1, 0.3, 0.5}, 
and /ir- = ±1 + O.li. 
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Figure 7: The angles 0i and 02 as functions of G (—10,10) and ej. G (0,10), for 9 G 
{30°, 150°}, /3 G {0.1, 0.3, 0.5}, and /i^ = ±1 + O.li. 
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Figure 7: Continued. 
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